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Abstract 

The electromagnetic El and Ml transitions in heavy quarkonia (cc,bb,cb) and the magnetic 
moment of the B^r are calculated within the framework of the covariant Blankenbecler- Sugar (BSLT) 
equation. The aim of this paper is to study the effects of two-quark exchange current operators 
which involve the QQ interaction, that arise in the BSLT (or Schrodinger) reduction of the Bethe- 
Salpeter equation. These are found to be small for El dominated transitions such as ip(nS) — > \ c j 7 
and T(nS) — * \bJ7t but significant for the Ml dominated ones. It is shown that a satisfactory 
description of the empirical data on El and Ml transitions in charmonium and bottomonium requires 
unapproximated treatment of the Dirac currents of the quarks. Finally, it is demonstrated that many 
of the transitions are sensitive to the form of the QQ wavefunctions, and thus require a realistic 
treatment of the large hyperfine splittings in the heavy quarkonium systems. 
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1 Introduction 



The radiative decays of heavy quarkonia (cc, bb, cb) have drawn much interest, as they can provide 
direct information on both the heavy quarkonium wavefunctions and the QQ interaction. Theoretical 
descriptions range from potential models - H, QCD sum rules || Heavy Quark Effective Theory 
(HQET) [^| to Non-Relativistic QCD (NRQCD) ||. As experimental data of a reasonable quality exists 
for a number of transitions in the cc and bb systems | fl0| [HJ , a fair assessment of the quality of model 
predictions is already possible. The measured 7 decays in the charmonium (cc) system include the 
El transitions Xcj — * J/i 3 ! and ip' ~ > XcJli as well as the spin-flip Ml transitions J /if) — > r) c j and 
i\) ' — > rjc'j. The situation in the bottomonium (bb) system is, however, less satisfactory as the total 
widths of the Xbj states are not known, and none of the spin-flip Ml decays observed. 

Previous calculations of the widths for El transitions in heavy quarkonia jl| |[ [l2| have acheived 
qualitative agreement with experiment. However, the situation concerning the Ml transitions has re- 
mained unsatisfactory for a long time fl3|| as the width for J/ip — > rj c j has typically been overpredicted 
by a factor of ~ 3. Calculations of Ml widths using the nonrelativistic Schrodinger equation |3|, |l4j have 
demonstrated that the Ml transitions in charmonium are sensitive both to the relativistic aspects of the 
spin-flip operator as well as the Lorentz structure of the QQ interaction. The results obtained in ref. |L4| 
suggest that a scalar confining interaction may explain the observed width of ~ 1 kcV for J/tp — > r] c j, 
provided that an unapproximated expression for the single quark spin-flip operator is used. However, 
as many of the Ml transitions are very sensitive to the exact form of the spin-flip operator and the QQ 
wavefunctions, a more realistic model for the wavefunctions of heavy quarkonia is called for. 

This paper reports a comprehensive calculation of the El and Ml transitions in heavy quarkonia, 
as well as the magnetic moment of the B^r (cb, be) meson, within the framework of the covariant 
Blankenbecler-Sugar (BSLT) equation Jl(j]. The QQ interaction Hamiltonian is formed of scalar con- 
fining and vector one-gluon exchange (OGE) components, in addition to a small instanton induced 
interaction |T^ ]. The hyperfine components of the QQ interaction have been fully taken into account, 
which is found to be important for several El and Ml transitions, in line with the conclusion reached in 
ref. [jl8| . On the other hand, the elimination of the negative energy components of the Bethe-Salpeter 
equation leads to the appearance of two-quark transition operators |]l9| . These exchange current opera- 
tors, which depend explicitly on the Lorentz coupling structure of the QQ interaction, have been shown 
to give large corrections to the single quark transition operators JlJ, . 

The aim of this paper is to evaluate the exchange current contributions to the electric and magnetic 
dipole operators for the QQ interaction Hamiltonian described above. The contributions from exchange 
charge operators to the El transition rates are shown to be highly suppressed by the large masses of the 
charm and bottom quarks. However, it is found that the exchange magnetic moment operator associated 
with the scalar confining component of the QQ interaction gives a large contribution, which makes it 
possible to bring the calculated width for J /tp — > 7^7 into agreement with the observed width of ~ 1 keV. 
It is also found that the empirical width of ~ 1 keV for the forbidden transition ip ' — > r/ c j may be similarly 
explained if a relativistic single-quark spin-flip operator is employed togehter with wavefunctions that 
model the spin-spin splitting of the ip states. This provides evidence in favor of a dominant scalar Lorentz 
structure for the effective linear confining interaction, as an effective vector interaction gives a vanishing 
contribution to the spin-flip operator for Ml transitions in the cc and bb systems [T^ |. On the other hand, 
for the B~- the OGE interaction contributes an exchange magnetic moment operator that counteracts 
that from the scalar confining interaction. 

The layout of this paper is as follows: Section 2 presents the transition operators for El and Ml 
decay, while section 3 deals with the Hamiltonian model and QQ wavefunctions, along with formulas for 
the El and Ml widths. Section 4 presents the numerical results for the radiative decays and the B^r 
magnetic moment, and section 5 contains a discussion of the obtained results. 
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2 The Electric Dipole and Magnetic Moment Operators 



2.1 The Charge Density and Electric Dipole Operators 

From the 5*-matrix element for one-photon emission by a two-quark system, 

Sfi = -i (f\ 8162 e M J M \i) , (1) 

where 8\ >2 are four-momentum conserving delta functions for quarks 1 and 2, and J M is the current 
operator of the two-quark system, the electromagnetic transition amplitude for QQ systems is obtained, 
in the impulse approximation, as 

T fi = - J d 3 nd 3 r 2 <p f (n, r 2 ) e ■ [e*** n(q) + e*** j 2 (q)] <Pi(n,r* 2 ), (2) 

where q and i denote the momentum and polarization of the emitted photon, respectively, while ipi 
and iff denote the orbital wavefunctions of the initial and final heavy quarkonium states. In the above 
equation, ji and j 2 denote the single quark current operators of quarks 1 and 2, respectively. Note that 
the current operator j(q) corresponds to the quantity in square brackets in eq. (|J). The amplitude for 
a 7 transition is then 



Tfi = i\q\ J d rid r 2 ip}(n, r 2 ) e ■ d (ri, r* 2 ) <Pi(fu ^2), (3) 
where the dipole operator d{fi,f 2 ) is of the form 

d(n,f 2 ) = [ d 3 r'e^ p 'r' p(f', r u r 2 ). (4) 



In order to distinguish the above model from the rigorous El approximation, it will be referred to 
as "dynamical" throughout this paper. The El approximation is then obtained in the limit q — > 0. 
If contributions to the charge operator p(f',r), that are proportional to higher powers of the photon 
momentum are to be included, then the usefulness of the dynamical model is apparent. More significantly, 
the dynamical model allows the recoil of the heavy meson to be accounted for. 

The charge density operator p(f') contains, in addition to the single quark contribution p sq , an 
exchange part p cx , which arises from processes that are illustrated by the diagrams shown in Fig. [I]. 
These contributions arise from the elimination of the negative energy components in the reduction of the 
Bethe-Salpetcr equation to a BSLT (or Schrodinger) equation. An obvious constraint is that two-quark 
contributions to the charge density should have vanishing volume integrals. 

Consider first the single quark charge operator p sq (f\ r) — pi(r', r\) + p 2 (r', f 2 ). The corresponding 
dipole operator may be expressed as 

d sq (ri, r 2 ) = Qi n e^rri + q 2 ^ e «>-ft ; (5) 

where Q\ and Q 2 are the appropriate electric charges of the constituent charm and bottom quarks. Note 
that Q\ is taken to be the charge of the heavy quark, while Q 2 denotes that of the heavy antiquark. 
Rclativistic modifications to the above expression will be considered further on. Insertion of eq. (|^) into 
ec l- (H|) yields the amplitude 

Tfi = i \q f \ (2tt) 3 5\Pf -Pi- q f ) f d 3 r <p* f {f) e ■ d(f) ^(r), (6) 
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where r — r\ — r-i, to which the contribution from the single quark dipole operator (||) is 



d sq {f) 



Qim 2 - Q2ITI1 



mi + 7712 



r e 



iqff/2 



(7) 



where r = fi —r%. Again, the El expression || is obtained by dropping the exponential in eq. (Q) and 
setting <ff = in the delta function in eq. (^). Note that for the cc and bb systems, the factor in brackets 
in eq. (J/]) reduces to the charge of the charm and bottom quark, respectively. 







Figure 1: Relativistic Born diagrams for photon emission by a heavy constituent quark. The lower 
diagrams describe the negative energy components of the upper diagrams, and can be obtained from 
the latter by separation of the intermediate quark propagators p a and pb into negative and positive 
energy components. Note that similar diagrams describe photon emission by the heavy antiquark. The 
exchange charge operators that correspond to the above diagrams have been calculated for different 
Lorentz coupling structures of the interaction V in ref . , which in the case of the QQ interaction will 
contain scalar confining and vector OGE components. 



Consider next the two-quark exchange charge operators from the Born diagrams given in Fig. ^. If 
those operators are decomposed according to p ex (r fi, fa) = p C xi(r' , r*i)+/0 G x2(r ', f 2 ), then the exchange 
charge contribution to the two-quark dipole operator 

4x (fi , f 2 ) = / d 3 r'e^ - F 'r' p cx (f ' , n , r 2 ) (8) 
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may be expressed as 
4a (r 1 )=f 1 e^' f?1 



d3k2 e-^-> exl (g>,fc 2 ) 



(2tt) 3 



lim 



rf 3 fc 2 

(2tt) s 



>exl (<f, fe) 



(9) 



from which the El approximation can again be obtained by setting qf — > 0. 

Having established eq. (H), the relativistic modifications to the single quark charge operator as well 
as the two-quark exchange charge operators from Fig. |l| may be considered. To second order in v/c, the 
single quark charge operator may be expressed in the form |2l|| 



Psq — Qi 



1L 

8m 2 



4m 2 



2), 



(10) 



where it is understood that the contribution from quark 2 is obtained by replacing the indices accordingly. 
The second term on the r.h.s. is the relativistic Darwin-Foldy term. It will be shown that the effect of this 
term is very small because of the large masses of the heavy constituent quarks. Note that this expansion 
is justified by the small coefficient of the q 2 term; It has been shown in ref. Hi that such an expansion 
cannot be used for the magnetic moment operator. The spin-orbit term in eq. (|l0|) is linear in the photon 
momentum q and will be left out in this work because of its smallness. For transitions between S-wave 
states, this term vanishes entirely Note that in the El approximation, the contribution from the 

Darwin-Foldy term to the dipole operator will likewise vanish. 

The exchange charge density operators that are associated with the QQ interaction have been cal- 
culated in ref. [ pl| , where the appropriate operators were extracted for different Lorentz invariants for 
systems composed of quarks with equal masses. When generalized to the case of unequal quark masses, 
the required operators are obtained as 



Pc 



Qi 

Amf 



q 2 V c (k 2 ) 



Q 2 
4m| 



q 2 V c {k 1 



for the scalar confining interaction, and 



Qi 

Am 2 



q- h 

mi 



2 q ■ k 2 Si ■ S 2 

3 m 2 



V g (k 2 



Q 2 

4m| 



q- ki 

m 2 



2 q- ki Si ■ S 2 

3 mi 



V g (ki 



(11) 



(12) 



for the vector coupled OGE interaction. In the above expressions, V c and V g denote the momentum space 
forms of the confining and OGE interactions, respectively. Note that the terms in eqs. ( |lT| ) and (|l^) that 
depend on k 2 correspond to the contribution p ex i in eq. (|9|) and vice versa. In the El approximation, 
the contribution to the dipole operator from eq. ( Jl~T| ) vanishes, while in the dynamical model, evaluation 
of eq. @ yields, analogously to eq. (^), 



< Conf 



(f) 



m 2 



Qi 

Amf mi + m 2 



Q2 



mi 



4m| mi 



m 2 



f V c (r) e lif f f/2 : 



(13) 



where the scalar confining interaction is of the form V c (r) — cr. Note that the second term in eq. ^ 
vanishes because q ■ i = 0. On the other hand, the OGE expression (12) gives a contribution, the 
dominant term of which is 



Qi 


( 1 ! ^Si-S 2 \ 




( 1 ! 1Si-S 2 


Ami 


\mi ' 3 m 2 J 


Am 2 


\m 2 3 mi 



dV g (r) 
r dr 



Jqff/2 



(14) 



Here V g (r) denotes the form of the OGE interaction in configuration space, and is here taken to be the 
Fourier transform of eq. (|32"|). This choice allows the inclusion of the running coupling of QCD. Note 
that eq. dlj) gives a non-vanishing contribution also in the El approximation. 
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2.2 The Current Density and Magnetic Moment Operators 

In the impulse approximation, the spin-flip magnetic moment operator for Ml transitions between S'-wave 
heavy quarkonium states may be obtained from the amplitude 



T 



fi 



{2n) i 5 i {Pf~P l -q f ) / d A ry* f {r)i 



,i(?-r/2 



-iq-r/2 



(15) 



where r — f\ — r*2- The corresponding matrix element for an Ml transition may be written in the form 



(16) 



M fi = i d d r ip* f (r) qxi-fl st ^(r ), 



where /2 s f denotes the spin-flip part of the standard magnetic moment operator. In addition to the single 
quark current operators in eq. (|l^), the negative energy Born diagrams in Fig. [^contribute a two-quark 
current operator J^ x . Decomposition of the single quark current into contributions from quark 1 and 2 
according to j sq — ]i + ]2 yields the magnetic moment operator 



lim 



-iq-r/2 



(17) 



in the nonrelativistic impulse approximation (NRIA). However, previous work has demonstrated that 
the static magnetic moment operators of the baryons are significantly modified by the canonical boosts 
of the constituent quark spinors (2^, |24|. In ref. |Q, it was shown that the static spin- flip magnetic 
moment operators for Ml decay of QQ states are also significantly affected, despite the large masses of 
the charm and bottom constituent quarks. 

The matrix element that corresponds to eq. ( |l6| ) in the relativistic impulse approximation (RIA) may 
be obtained as 



M*f l = i 



d 3 P 
(2^)3 



d 3 rd 3 r' e lp - {F '~ ?) <p* f (f) qx e-flf^(P) Vi {r), 



(18) 



where the final and initial state coordinates r ' and r are defined as r\ — f' 2 and r\ — f*2 respectively. 
In eq. (f[8|), the momentum variable P is defined as P = (p' + p)/2, where p' and p are the relative 



momenta in the representation p\ — Pi/2 + p, p 2 — Pi/2 — p and p[ = Pf/2 + p' , p' 2 = Pjj2 — p'. 



The relativistic single quark magnetic moment operator that appears in the matrix element (|18|) can be 
obtained from eq. (|17| ) by the substitution r — > (r ' + r)/2. 

In the nonrelativistic case, the spin-dependent part of the single quark current operator J^ q = Ji + J2 
is of the form 



J 



— {&x +a 2 ) x q 



Qi , Q2 



2mi 2m 2 



— ((?! -a 2 ) X q 



Qi Q2 



2mi 2m 2 



(19) 



The first term, which vanishes for equal mass quarkonia, describes the magnetic moment of the two-quark 
system whereas the second term is the spin-flip operator for Ml decay in the nonrelativistic impulse 
approximation (NRIA). The corresponding spin-flip operator in the relativistic impulse approximation 
(RIA) has been calculated in refs. J|4|, [2(i|| , and may for transitions between S'-wave states be expressed 
as 



P sq 



2mi 



-a 



2m 2 



2 , 7 

J2 



(<?1 - CT 2 ), 



si _ 

Jl 3Ei 



2+ K 



(20) 



with Ei = ^fW 



It is apparent from the above expression that the relativistic treatment will 



effectively weaken the NRIA result. 
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The relativistic Born diagrams in Fig. [j] contribute significant two-quark exchange currents that 
give rise to two-quark magnetic moment operators |23]. This situation is akin to that for the mag- 



netic moments of the baryons 123L 24 , but in that case additional complications arise from flavor de- 



pendent meson exchange interactions. Decomposition of the exchange current operator according to 
Jcx(<f, k\,k 2 ) — ]cxi(q, k 2 ) + jcx2(<T, ki) yields the exchange magnetic moment operator 



Me 



lim 

9-0 



V x 



Jq-f/2 



d 3 k : 

W) 



2 e- lk ^j cxl (q,k 2 ) 



-iq-r/2 



(21) 



However, the exchange magnetic moment operators turn out to be difficult to calculate directly from the 
above equation. Instead, they can be conveniently extracted from the expression j22| 



Me 



lim 

?-»0 



d A k 



e-^Va x 



Jcxl 



■ Jex2 



2 fc 



(22) 



where k 2 — k + q/2 and k\ = —k + q/2. By means of eq. (p2[), it is now possible to consider the two-quark 
current operators for the scalar confining and vector OGE interactions, as calculated from the diagrams 
in Fig. [l]by refs. [2^]. The two-quark current operator associated with the scalar confining interaction 
is of the form 



Q\Pi , QIP2 



+ -((?!+ 0*2) x q 



2m{ 



Q*2 

2ml 



+ 0(^1 - ^2) x q 



2m{ 



2ml 



(23) 



in which case the two-quark magnetic moment operator is most conveniently computed using eq. 
as the spin part of eq. (^3|) depends explicitly on the photon momentum q. In the above equation, the 
variables Q\ and Q\ are defined as Q\ = V c (k2)Qi and Q\ — V c {k±)Q2, respectively. The corresponding 
current operator for the OGE interaction may be expressed as 



Qi 



ia\ x k 2 , 2P 2 + iS 2 x k 2 



2m\ 



2m\m 2 



+ Q2 



ia 2 x ki 2P\ + io\ x ki 



2ml 



2m\m 2 



(24) 



with Q\ = V g (k 2 )Qi and Q 2 — V g (ki)Q 2 . As the above equation depends only on fci and k 2 , the OGE 
magnetic moment operator is most conveniently calculated using eq. ([22]). By Fourier transformation, 
the resulting magnetic moment operators for transitions between S'-wave quarkonium states may be 
obtained as [EQl 



Me C x° nf 



02 



((?i - a 2 ) + 



Qi , Q2 



for the scalar confining interaction, and 



Me 



eV g (r) 



Q2 

™2 



m\m 2 



(<?i - a 2 ) + 



1)17 



III') 



((71 + CT 2 ) 



02 

™2 



Q1+Q2 



m\m 2 



((?! + (7 2 



(25) 



(26) 



for the OGE interaction. For equal constituent quark masses, eqs. ([25|) and ([26|) reduce to the expressions 
given in ref. p2| . Note that the presence of a spin-flip term in the OGE operator ( |2rj ) is solely a 
consequence of the difference in mass between the constituent quarks, and will thus not contribute to the 
Ml decay widths of the charmonium and bottomonium states. Similarly, the terms that are symmetric 
in the quark spins vanish for equal mass quarkonia. However, in the case of the system, these terms 
will contribute to the magnetic moment of the cb system. Also the spin-flip Ml decays in the B^r system 
will receive a contribution from the OGE operator. 
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3 Wavefunctions and Decay Widths 



3.1 Hamiltonian model 

The interaction Hamiltonian employed in this paper for the QQ interaction contains contributions from 
the scalar confining, the vector one-gluon exchange (OGE), and the instanton induced interaction of 
ref. 111711 • This Hamiltonian is thus of the form 



-Hint = Vbonf + Vfoge + V[nst , 



(27) 



and has been employed together with the covariant Blankenbecler-Sugar (BSLT) equation. That equation 
may be expressed as an eigenvalue equation of the form 



{H + fli„ t )*Blm(r ) = e(E, Mq, M Q )^ nlm (r ), 



(28) 



where Ho denotes the kinetic energy operator of the nonrelativistic Schrodinger equation and the eigen- 
value e(E, Mq, Mq) is a quadratic mass operator Jl5| which is related to the energy E of the heavy 
quarkonium state according to 



[E 2 - (M Q + Mq) 2 } [g - (Mq Mq) 2 } 
8fiE 2 



(29) 



where fi is the reduced mass of the quark-antiquark system. As the details of the interaction Hamil- 
tonian (27) have already been described in ref. only the main points will be repeated here. The 
interaction operators V in eq. (p7|), which in general are nonlocal, may be obtained from the QQ irre- 
ducible quasipotential V according to 



V(p',p) = 



Q 



Mr 



W(p') 



V(p',p) 



W(p) 



(30) 



where the function W is defined as W(p) — Eq(p) + Eq(p) with Ei(p) = sjM'l +p 2 . Note that in 
the Born approximation the quasipotential V is set equal to the QQ invariant scattering amplitude T, 
whereby a constructive relation to field theory obtains. 

The OGE interaction in momentum space can be parameterized [p5| in terms of the strong coupling 
a s (k 2 ) according to 



^Oge(fc) 



167T a s (k 2 ) 
~3 



as (p) = ™L l n -i 
' 27 



k 2 + Am 2 
~1\ 2 

J1 "QCD 



(31) 



Here Aqcd denotes the QCD scale of ~ 250 MeV, and rru is a dynamical gluon mass which determines 
the low-momentum behavior of a s . Application of eq. ([nj) together with eq. (|3^) yields the central and 
spin-dependent potential components of the OGE interaction, and are given in ref. fl5|| . For example, 
if all higher order nonlocalities are dropped, the central Coulomb component of the OGE interaction is 
modified to 

4 2 f°° Mr, Mr, ( Mq + M n \ 

Vo ge (r) = ~- dfcjofcr— — \ Q )a s (k 2 ), (32) 

where the factors eQ and eg are defined as eQ = \J Mq + k 2 /4 and eg = ^Mq + ^ 2 /4. If a s is taken 
to be constant, then the above form reduces to the Coulombic potential suggested by pcrturbative QCD 
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in the limit Mq — * oo. On the other hand, the Fourier transform of a linear confining potential of the 
form Vconf = cr may be expressed as 



Vr (D=lim l-^Le-^ 8 ™ 



a^o J (2tt) 3 fc 4 



4A 2 fc 4 k 4 



(A 2 + fc 2 ) 3 (A 2 + k 2 ) 2 



(33) 



For the purpose of calculating the spin-orbit term associated with the scalar confining interaction, the 
limit A — > may be taken directly to yield Vconf (k ) — Sirc/k 4 . Finally the instanton induced interaction 
is expressed as 

AM Q AM Q r~ { M Q + M Q \ MqMq 

V[nst (r = / dk k j (kr) ■ ^ (34) 

4 « Jo \ e Q + e Q J e Q e Q 

where the notation is similar to that employed in ref . |0] . The parameter AMq denotes the mass shift 
of the heavy constituent quark due to the instanton induced interaction, which for a charm quark is of 
the order ~ 100 MeV [ ft7| . The parameter n represents the instanton density, which is usually taken to 
be ~ 1 fm -4 . Note that in the limit of infinitely heavy constituent quarks, eq. (Q) reduces to a delta 
function. The above smeared out form is convenient since it allows for direct numerical treatment of the 
instanton induced interaction with the differential equation (pgj). 

The wavefunctions needed for the calculation of the El and Ml widths of the heavy quarkonia are 
thus taken to be the solutions to eq. (^) obtained in ref. As the spin-spin interaction in the 

charmonium system is strong enough to produce a J/tp — r/ c splitting of ~ 120 MeV, then the respective 
radial wavefunctions are likely to show marked differences. Thus the employment of spin-averaged 
wavefunctions and a perturbative treatment of the spin-dependent hyperfine interaction is undesirable. 
With this in mind, the hyperfine components of the QQ interaction have been taken fully into account, 
which is possible since they are nonsingular in the BSLT framework. In particular, it is expected that such 
effects should be significant for transitions which involve a change in the principal quantum number of 
the QQ state. The effect of the 5-wave spin-spin interaction on the radial wavefunctions is demonstrated 
by Fig. ^. The calculated spectra of the heavy quarkonia and the model parameters are given in Tables@ 
and ||, respectively. 



0.01 




Figure 2: The reduced radial wave- 
functions for the charmonium IS* 
and 2S states, from ref. (l5j. The 
differences between the spin singlet 
rj c and spin triplet tp wavefunctions 
are due to the short-ranged OGE 
spin-spin interaction. Note that the 
r axis has been made logarithmic in 
order to emphasize the short range 
part of the wavefunctions. 
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Table 1: Calculated and experimen- 
tal charnionium, bottonionium and 
B^r states rounded to the nearest 
MeV. The states are classified ac- 
cording to excitation number n, to- 
tal spin S, total orbital angular mo- 
mentum L and total angular mo- 
mentum J. Note that experimen- 
tal uncertainties are indicated only 
where they are appreciable. For a 
graphical plot of the cc and bb data, 
see ref. The cc and bb values 
are from ref. p"5| . The experimen- 
tal states correspond to the values 
reported by ref. fnl , except for the 
recently observed ||7) 7? c (25). The 
measured mass of the B^r was re- 
ported in ref. @ as 6.40±0.39 GeV, 
which is about ~ 100 MeV higher 
than the predicted 6308 MeV, and 
most other models give even lower 
masses for the B^r ground state ||. 
The predicted B^r spectrum agrees, 
however, very well with the QCD- 
inspired model of ref. (^9). Within 
the framework of QCD sum rules, 
the mass of the Bf has been pre- 
dicted to be ~ 6.35 GeV @. 





Ref. || 


Other sources 


M h 
M c 

Aqcd 
m g 
c 


4885 MeV 
1500 MeV 

260 MeV 
290 MeV 
890 MeV/fm 


4870 MeV @ 
1530 MeV |g] 

200-300 MeV j25[ 
m g > Aqcd pff 
912 MeV/fm Q 


(AM c f 
in 

(AM„) 2 
4n 

AM b AM a 
4n 


0.084 fm 2 
0.004 fm 2 
0.018 fm 2 


- 0.05 fm 2 @ 



Table 2: Constituent quark masses and parame- 
ters for the QQ interaction that have been used in 
the calculation of the spectra presented in Table [I]. 
The heavy masses are close to those preferred by 
ref. fl30T| , and in general in agreement with the val- 
ues of previous work. The form of the running 
QCD coupling a s is in line with the criteria of 
refs. 1 25, and the confining string tension c 
agrees well with previous calculations |?0[ [3l| . The 
strength of the instanton induced interaction for cc 
is comparable to the estimate given by ref. |l7|. 
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3.2 Widths for radiative decay 

The decay width for an El dominated transition of the type Xcj ~ * J 7^7 or "0' — > XcJl is given by 



S 



Si 



2J f + l 



M 



q a 



Mi 



\M Q \ 



\M 2 \ 2 



(35) 



where Jf is the total angular momentum of the final quarkonium state, and q is the momentum of the 
emitted photon. One would thus expect that the widths for ip' — > \c,J 7 with J = 0, 1, 2 would scale as 
1:3:5 respectively. In practice, however, this result is modified by the large hyperfine splittings in 
the L = 1 heavy quarkonia. The factor Sfi is defined as in ref. [[| and assumes the values Sfi — 1 for a 
triplet-triplet transition and Sfi — 3 for a singlet-singlet transition of the type h c — > r\d- On the other 
hand, the widths for transitions between P- and P-wave states are calculated according to 



45 



2J/ + 1 3 M f ... . 2 



5 



fi 



18 



i 



(36) 



where Jd and J p are the total angular momenta of the D- and P-wave states, respectively. Note that the 
triangularity of the 6-j symbol requires that | Jd — J p \ = 1 or 0. Consequently, transitions that change the 
value of J by more than one unit are forbidden. In eqs. (|3^) and (|36|), Mq and M2 denote radial matrix 
elements for S- and P-wave photon emission, respectively. The radial matrix element for 5-wave decay 
receives contributions not only from the impulse approximation, eq. (|?]), but also from the confinement 
and OGE operators (Il3|) and (fl4|). That matrix element may thus be expressed as 



A4q = / drruj 
'0 



(r)ui(r) j (y) 



(Q) IA + q 2 V c (r) (Q) t 



dV g (r) 
r dr 



(Q)< 



(37) 



where Ui and Uf are the reduced radial wavefunctions for the initial and final heavy quarkonium states. 
Similarly, the matrix element for P-wave decay is of the form 



poo 

M 2 = (Q)i D / drru f (r)ui(r) j 2 
Jo 



(38) 



The contribution from this matrix element is generally very small, and has therefore not been included 
in eq. fl36|). 

The impulse approximation charge factor (Q) IA , and the exchange charge factors (Q) c for the scalar 
confining interaction and (Q) g for the OGE interaction that appear in eqs. ( |37| ) and ( |38| ) are of the form 



(Q) 



IA 



Qi 1 



m 2 



8m^ J mj + rri2 



Q2 1 



3 



mi 



77T-2 / JTll + ?7l2 



(39) 



for the impulse approximation, where the quark charge operators have been multiplied with the Darwin- 
Foldy terms from eq. (|l(]), and 



(Q) c 



m 2 



Q2 



mi 



Am^ mi + ni2 Am^ mi + 1712 

Qi_ (J_ + 2 (Sf\a 1 -a 2 \S i ) \ _ J h _ ( L 2 (.S>|^-^ 7,) 
4m^ 3 



m 2 



Ami \TO2 



/7?1 



(40) 
(41) 



for the confinement and OGE exchange charge contributions, respectively. In the spin dependent terms 
of eq. (|4l|), Si and Sf denote the total spins of the initial and final quarkonium states. For triplet-triplet 
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and singlet-singlet transitions, (S/|ov <?2|5j) = +1 an d —3, respectively. The charge factor (Q) m that 
appears in eq. ( p^ ) is defined according to (Q) ID = lim^o (Q)ia- This is permissible since the Darwin- 
Foldy and exchange charge terms are very small compared to the dominant dipole contribution, which in 
itself is already insignificant because of the suppression by the j% function in the matrix element. There 
is thus no need to include these terms in the matrix element M2 for D-wave decay. 

In the El approximation, the recoil factor Mf/Mi vanishes by eq. (||), since in that case Pf = Pi. 
The decay width expression (J3q) thus reduces to 



r 



El 



2Jt + l , 
= 5/i— £—4 Co* limMc 
27 g->0 



(42) 



in the El approximation, which is similar to the expression given in ref. j5[. Here the "recoilless" 
q- value is given by q nr = Mi — Mf. Note that the OGE exchange charge contribution survives in the El 
approximation, whereas the contribution from the scalar confining interaction as well as the Darwin-Foldy 
terms are eliminated. 

The expression for the width of a spin-flip Ml transition between 5*-wave heavy quarkonium states 
can be written in the form 

Tmi = J^l^^lM.f, (43) 



2Si 



where M. 1 denotes the radial matrix element for Ml decay and Si is the total spin of the initial state. 
The radial matrix element for Ml decay consists of relativistic impulse approximation, scalar confining 
and OGE components, according to 



M- 



M 



RIA 



M 



Conf 



x° sc , 



(44) 



where the different matrix elements can be obtained from the corresponding spin-flip magnetic moment 
operators by dropping the charge e and the spinors u\ — 8%. From the relativistic matrix element (18) 
and the relativistic spin- flip magnetic moment operator (^of), one obtains 



M 



RIA 



2 POO roo />oc 1 

- / dr'r'utir 1 ) \ drr Ul (r) / dP P 2 - 
k Jo Jo Jo 4 



Qi 

mi ' 



fl 



9i 

rri2" 



& 



jo(r'P)jo(rP), (45) 



where the factors / z 7 are given by eq. ( p0| 
vector OGE interactions are of the form 



The matrix elements associated with the scalar confining and 



M 



Conf 

7 



druf(r) Ui(r) 



V c (r) 



Qi O2 



(46) 



for the scalar confining interaction, and 



M° sc = 



druf{r) Ui{r) 



V g (r) 



Qi Qi Qi — Q2 



TO1TO2 



(47) 



for the OGE interaction. 

It should be noted at this point that the appearance of two-quark matrix elements as given above 
is entirely due to the elimination in the impulse approximation of the negative energy components in 
the reduction of the Bethe-Salpeter equation to a Blankcnbecler-Sugar equation. These components do 
however contribute as transition matrix elements, and have to be included in the Blankenbecler-Sugar 
(or Schrodinger) framework as explicit two-quark current operators Jl9| , as illustrated in Fig. [l]. In 
e.g. the alternate Gross type reduction of the Bethe-Salpeter equation, these two-quark operators are 
automatically taken into account by the single quark transition operators. 
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3.3 Magnetic moment of the Bf 

Of the heavy quarkonium systems, only the spin-triplet Bf (cb, cb) states have a magnetic moment. That 
may be calculated from the spin-symmetric part of eq. When expressed in terms of the nuclear 

magneton /ijy, the magnetic moment of the B* + (cb) state is obtained as 



M c5 - — 



2 



1 



Mat, 



(48) 



where m p is the proton mass. The above expression is valid in the nonrelativistic limit. Since it is known 
that the magnetic moments of the baryons receive significant corrections from relativistic effects p3f , 
then a relativistic version of eq. ( |4§| ) is called for. In the relativistic impulse approximation, the magnetic 
moment of an S'-wave cb state is of the form 



dr dr' r r u(r) u(r') 



dPP 2 ^ 
3 



irn + —n 



i 



Jo (r'P)jo (rP) m, 



(49) 



where u(r) is the reduced radial wavefunction of the cb state. In eq. ([49|), the factors are defined in 
eq. (p0|). In addition to the relativistic impulse approximation, the exchange magnetic moment operators 
associated with the scalar confining and vector OGE interactions will also contribute to the magnetic 
moment of the B^ ■ The contribution from the scalar confining interaction is of the form 



Conf 



dru 2 (r) V c (r) 



1 



(50) 



where the confining interaction V c (r) is of the form V c (r) 
expressed as 



cr. Similarly, the OGE contribution can be 



^cb 



dr u 2 (r) V g (r) 



6 



m c m b 



MAT) 



(51) 



where the OGE potential V g (r) is given by the Fourier transform of eq. (J3T]) . The total magnetic moment 
of the B*^ is thus given by the sum of the RIA contribution (^) and the exchange contributions (50) 
and @. 



4 Numerical Results 

This section presents the numerical values of the widths and matrix elements for each El and Ml 



transition, as obtained using the model for the QQ spectra presented in section 3.1. Unless otherwise 
indicated, the El widths have been calculated using eqs. (|35|), (|36|), ( |37| ) and (|38|), while the widths for 
Ml decay correspond to eqs. ( fl3[ ) and (|44|). When the mass of one of the quarkonium states is not known 
empirically, then the splittings rather than the absolute values from Table [j] are used to determine the 
photon momentum q y . 

The widths for Ml transitions between S'-wave states in charmonium and bottomonium are given 
in Table |^, along with the associated impulse approximation and exchange current matrix elements. 
For cc and bb, the width for Ml decay is given by the sum of the impulse approximation and scalar 
confinement terms. In the case of the B^, the widths for Ml decay also receive a contribution from the 
OGE exchange current. The Ml widths and magnetic moments of the B^ 1 states are given in Tables [lC] 
and [l^, respectively. The widths of the El dominated transitions between low-lying cc and bb states are 
given in Tables |J, |[ ^, @ and [|. They have been calculated both for the dynamical model presented in 
section 2A and the rigorous El approximation. The computed widths for El dominated transitions in 
the Bf system are given, along with the matrix elements, in Tables ^, |ll| and [l2]. 
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Table 3: The Ml transitions between low-lying 5-wave states in the charmonium (cc) and bottomo- 
nium (bb) systems. Experimental data jl(J is available only for the J/ip — ► rj c -f and ifi' — > r/ c j transi- 
tions. Note that the empirical value for ip ' — > r\ c ^ is uncertain since the total width of the if) ' is poorly 
known. The quoted photon momenta q 1 have been obtained by combination of the empirical masses of 
the spin triplet states jnj with the splittings given by the Hamiltonian model of ref. Jll| in Table [l]. The 
Ml decays of the ^(35) state have not been included since the 35 states in charmonium lie above the 
threshold for DD fragmentation. 



Transition 






Matrix element [fm] 








Width 




NRIA 




RIA 


Conf 




NRIA 


RIA 


RIA+Conf 


J/ip-> Vcl 
<j 7 : 116 MeV 


4.356 


10 


-2 


3.762 • 10~ 2 


-8.724 


10" 


-3 


2.85 


2.12 


1.25 keV 

exp: 1.14 ±0.39 


1p ' -> Vcl 
q-y : 639 MeV 


3.985 


10 


-3 


-5.14 • 10~ 4 


2.826 


1(T 


-3 


3.35 


0.06 


1.13 keV 

exp: 0.84 ±0.24 


V> ' -> Vcl 
<j 7 : 46 MeV 


4.344 


10 


-2 


3.735 • 10~ 2 


-1.870 


io- 


-2 


0.18 


0.13 


0.03 keV 


v'c -> J/i^i 

q-y : 502 MeV 


-4.271 


10 


-3 


-7.584- 10" 3 


5.206 


10" 


-3 


5.89 


18.6 


1.83 keV 


T — ► ?/f,7 
<j 7 : 59 MeV 


-6.71 


10 


-3 


-6.39-10- 3 


2.41 


10" 


-4 


9.2 


8.3 


7.7 eV 


T(25) - Vbl 

q-t : 603 MeV 


-3.94 


10 


-4 


-1.44- 10~ 4 


-8.76 


10" 


-5 


31.8 


4.3 


11.0 eV 


T(25) ^7/ b (25) 7 
<? 7 : 25 MeV 


-6.70 


10 


-3 


-6.39- 10~ 3 


5.45 


10" 


-4 


0.70 


0.64 


0.53 eV 


r) b (2S) -» T 7 

q 7 : 530 MeV 


4.18 


10 


-4 


6.30- 10~ 4 


-1.31 


10" 


-4 


71.5 


162 


102 eV 


T(35) ^r/ b (35) 7 

q-y : 16 MeV 


-6.70 


10 


-3 


-6.35 • 10~ 3 


8.02 


10" 


-4 


0.18 


0.16 


0.13 eV 


T(35) ^r/ b (25) 7 

<j 7 : 350 MeV 


-3.59 


10 


-4 


-1.11 • 10~ 4 


-1.55 


10" 


-4 


5.3 


0.5 


2.9 eV 


T(35) -> r/ b 7 

<j 7 : 910 MeV 


-2.10 


10 


-4 


-6.59- 10" 5 


-3.77 


io- 


-S 


30.2 


3.0 


7.3 eV 


77 b (35)^T(25) 7 

q~, : 311 MeV 


3.96 


10 


-4 


6.05 • 10~ 4 


-2.25 


io- 


-4 


13.7 


32.0 


12.6 eV 


776(35) - T 7 

<j 7 : 842 MeV 


2.05 


10 


-4 


3.02 • 10~ 4 
14 


-4.68 


io- 


-S 


68.7 


149 


106 eV 



Table 4: The El dominated transitions between low-lying states in the charmonium (cc) system, together 
with the empirical data given by ref. |jfj]| . The column "IA" contains the matrix element ( |37| ) in the 
impulse (single quark) approximation, while in the column labeled "DYN", the exchange charge contri- 
butions have been included. The columns "El" contain the matrix element and the 7 width in the El 
approximation. The g 7 values, as given above, have been rounded to the nearest MeV, and correspond 
wherever possible to the empirical data in Table 0. Further El transitions can be found in Table @. 



Transition 




M [fm] 




M 2 


fm] 




Width 




IA 


DYN 


El 








El 


DYN 


Xc2 -> J/ipm 

q-y : 429 MeV 


0.2389 


0.2442 


0.2632 


7.145 


• 10 


-3 


558 keV 


343 keV 

exp: 389 ±60 


Xcl -> J/ipH 

<? 7 : 390 MeV 


0.2464 


0.2519 


0.2673 


5.729 


• 10 


-3 


422 keV 


276 keV 

exp: 290 ± 60 


XcO — * J 1^1 
<j 7 : 303 MeV 


U.ZOOD 


U.Z01Z 


n 07m 


o.o4o 


1 


-3 


196 keV 


144 keV 

exp: 165 ± 40 


'0 ' —* XcO 7 

q-y : 261 MeV 


— U.ZDoO 


— U.ZDoD 


n os/in 

— U.Z04U 


— D.1UD 


1 n 


-3 


44.6 keV 


33.1 keV 

exp: 26.1 ± 4.5 


V>' — > Xci 7 

q-y : 171 MeV 


n qi OA 
— U.olZD 


— U.ol/D 


O QOOO 


Q HOC 


1 


-3 


45.8 keV 


38.7 keV 

exp: 25.2 ± 4.5 


ip ' — *• Xc2 7 
<j 7 : 127 MeV 


Q/i/in 

— U.G44U 


Q/i/io 


O Q/1 so 

— u.o4oy 


1 Q71 
— l.Oi 1 


1 

• lu 


-3 


37.1 keV 


33.1 keV 

exp: 20.4 ± 4.0 


h c —> rj c 7 

<j 7 : 493 MeV Sfi : 3 


n oooc 


u.zuy 1 


n ooco 


7 Q77 
1.0/ i 


1 
• lu 


-3 


661 keV 


370 keV 


<j 7 : 125 MeV S fi : 3 


— 0.3420 


n /i a 
— 0.3424 


— 0.3465 


— l.olo 


■ 10 


-3 


61.5 keV 


55.0 keV 


^(35) -> XcO 7 

? 7 : 577 MeV 


-0.0456 


-0.0450 


-0.0199 


0.926 


• 10 


-2 


2.69 keV 


9.86 keV 


^(35) -> Xci 7 

<j 7 : 494 MeV 


-0.0306 


-0.0298 


-0.0033 


1.016 


• 10 


-2 


0.13 keV 


9.57 keV 


V'(3'S') -> Xc2 7 

<j 7 : 455 MeV 


-0.0168 


-0.0161 


0.0123 


1.099 


• 10 


-2 


2.38 keV 


5.75 keV 


V(35)^Xco(2P)7 

q-y : 153 MeV 


-0.4315 


-0.4315 


-0.4497 


-7.344 


• 10 


-3 


21.3 keV 


17.8 keV 


V(3S)^Xcl(2P) 7 
q-y : 125 MeV 


-0.4860 


-0.4861 


-0.4995 


-5.399 


• 10 


-3 


42.6 keV 


37.3 keV 


V>(3S)^Xc 2 (2P) 7 
q~, : 109 MeV 


-0.5280 


-0.5283 


-0.5391 


-4.367 


• 10 


-3 


53.7 keV 


48.2 keV 
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Table 5: The El dominated transitions between low-lying states in the bottomonium (bb) system, to- 
gether with the empirical data given by ref. j[o|. The column "IA" contains the matrix element (37) 
in the impulse (single quark) approximation, while in the column labeled " DYN" , the exchange charge 
contributions have been included. The columns "El" contain the matrix element and the 7 width in the 
El approximation. The photon momenta q 1 have been obtained as for Table || Further El transitions 
that involve the Xbj (2P) states can be found in Table 13. 



Transition 




Mo [fm] 




M-2 


[fm] 




Width 




IA 


DYN 


El 








El 


DYN 


Xb2 -> T 7 

q 1 : 443 MeV 


-0.0721 


-0.0723 


-0.0743 


-7.76 


• 10 


-4 


42.7 keV 


36.0 keV 

Br: 22 ±3% 


Xbi -> T7 

q 1 : 424 MeV 


-0.0731 


-0.0733 


-0.0751 


-6.96 


• 10 


-4 


38.1 keV 


32.5 keV 

Br: 35 ± 8 % 


Xbo -» T 7 

<? 7 : 392 MeV 


—0.0742 


—0.0743 


A ATC A 

—0.0759 


C OA 

—5.80 


1 A- 

■ 10 


-4 


30.7 keV 


26.6 keV 

Br: < 6% 


'V /o n\ 

T(25) -> Xbo7 

<? 7 : 162 MeV 


0.0935 


0.0935 


A AA A A 

0.0942 


2.90 


• 10 


-4 


1.07 keV 


1.01 keV 

exp: 1.7 ±0.5 


'V /o n\ 

<? 7 : 129 MeV 


0.1007 


0.1007 


0.1012 


2.01 


• 10 


-4 


1.88 keV 


1.80 keV 

exp: 3.0 ± 0.7 


'V /o o\ 

<? 7 : 109 MeV 


0.1063 


0.1063 


A 1 A/"*T 

0.1067 


1.53 


1 A 

• 10 


-4 


2.10 keV 


2.03 keV 

exp: 3.1 ± 0.7 


<j 7 : 485 MeV S/j : 3 


— 0.0do4 


— 0.0od4 


—0.0685 


T AO 

— 7.98 


1 A- 

■ 10 


-4 


47.9 keV 


39.7 keV 


776(25") ->• /l&7 

<j 7 : 100 MeV S/j : 3 


0.1066 


0.1066 


0.1069 


1.20 


• 10 


-4 


2.86 keV 


2.77 keV 


T(3S) Xbol 

q 1 : 483 MeV 


0.0068 


0.0068 


0.0039 


-1.13 


• 10 


-3 


0.05 keV 


0.15 keV 


T (3S') X617 

g 7 : 452 MeV 


0.0036 


0.0036 


0.0005 


-1.20 


• 10 


-3 


2.2 eV 


0.11 keV 


T(3S) -> x 6 2 7 

g 7 : 433 MeV 


0.0007 


0.0007 


-0.0025 


-1.26 


• 10 


-3 


0.08 keV 


0.04 keV 


T(35)^x 6 o(2P)7 

<? 7 : 122 MeV 


0.1505 


0.1504 


0.1520 


6.18 


• 10 


-4 


1.19 keV 


1.14 keV 

exp: 1.4 ±0.3 


T(35)^xw(2P)7 

<j 7 : 100 MeV 


0.1614 


0.1613 


0.1624 


4.36 


• 10 


-4 


2.20 keV 


2.12 keV 

exp: 3.0 ±0.5 


T(35)-^x 62 (2P)7 

q-y : 86 MeV 


0.1699 


0.1699 


0.1707 


3.38 


• 10 


-4 


2.57 keV 


2.50 keV 

exp: 3.0 ±0.6 
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Table 6: The El transitions from the lightest spin-triplet D-wave states in charmonium (cc). The 
labeling of the columns is as for Tables and and the photon momenta q 1 have been calculated from 
the D- wave masses predicted in Table [jj The statistical factors Sfi are given by eq. (|3^). Note that a 
good experimental candidate for the 3 D\ state is the ip (3770) resonance. 



Transition 






A A \C 1 

M [fmj 




Width 






IA 


DYN 


El 


El 


DYN 


3 D 3 -V X c2 7 

q„, : 227 MeV S fi : 


18/25 


0.4164 


0.4194 


0.4353 


243 keV 


192 keV 


3 D 2 -> Xc2l 

q-f : 221 MeV Sfi : 


9/50 


0.4188 


0.4219 


0.4367 


56.5 keV 


45.2 keV 


3 D 2 -> Xci 7 

9 7 : 263 MeV Sfi : 


9/10 


0.3920 


0.3953 


0.4145 


262 keV 


198 keV 


3 £>i -> Xc2 7 

<j 7 : 206 MeV Sfi : 


1/50 


0.4216 


0.4246 


0.4372 


5.06 keV 


4.13 keV 


3 Dt -> Xci7 

<j 7 : 248 MeV <S/ 4 : 


1/2 


0.3963 


0.3997 


0.4164 


123 keV 


94.9 keV 


3 L>1 -> Xc0 7 

g 7 : 336 MeV Sfi : 


2 


0.3578 


0.3619 


0.3889 


370 keV 


251 keV 



Table 7: The El transitions from the lightest spin-triplet D-w&ve states in bottomonium (bb). The 
photon momenta q 1 and statistical factors Sfi have been calculated as for Table ||. 



Transition 








Mo [fin] 






Width 








IA 


DYN 


El 


El 


DYN 


3 D 3 -» Xb2 7 

q 1 : 242 MeV Sfi : 


18/25 


-0 


1270 


-0.1271 


-0 


1291 


24.5 keV 


22.3 


keV 


3 D 2 -> Xb2 7 

q 1 : 233 MeV Sfi : 


9/50 


-0 


1275 


-0.1276 


-0 


1295 


5.51 keV 


5.04 


keV 


3 D 2 -» Xbi 7 

q-y : 253 MeV S fi : 


9/10 


-0 


1228 


-0.1229 


-0 


1250 


19.6 keV 


17.8 


keV 


3 D\ -> Xb2l 

g 7 : 223 MeV Sfi : 


1/50 


-0 


1280 


-0.1280 


-0 


1297 


0.54 keV 


0.50 


keV 


3 Di -> Xbil 

q-y : 243 MeV Sfi : 


1/2 


-0 


1234 


-0.1235 


-0 


1254 


9.75 keV 


8.89 


keV 


3 Di -> XbOl 

q-y : 275 MeV Sfi : 


2 


-0 


1173 


-0.1174 


-0 


1196 


17.3 keV 


15.5 


keV 
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Table 8: The El transitions from the 2P states in charmonium (cc) and bottomonium (bb). The labeling 
of the states and the photon momenta g 7 for the transitions are as for Tables | and §. 



Transition 




Mo [fm] 




M-2 


[fm] 




Width 




IA 


DYN 


El 








El 


DYN 


Xc2{2P) -> tp 7 

<j 7 : 236 MeV 


0.3702 


0.3735 


0.4007 


i no 

1.09 


■ 10 


-2 


194 keV 


144 keV 


/o 7~>\ / / 

Xci(2P) -» V 7 

q 1 : 220 MeV 


0.3996 


0.4031 


0.4271 


9.60 


• 10 


-3 


178 keV 


137 keV 


<j 7 : 193 MeV 


0.4321 


0.4357 


0.4543 


7.43 


• 10 


-3 


133 keV 


108 keV 


Xc2 (2P)^ J/^7 

<j 7 : 745 MeV 


0.0660 


0.0693 


0.0526 


-6.49 


• 10 


-3 


133 keV 


132 keV 


Xci(2P)^ J/^7 

<j 7 : 731 MeV 


0.0549 


0.0581 


0.0380 


-7.83 


• 10 


-3 


65.3 keV 


90.1 keV 


g 7 : 707 MeV 


0.0387 


0.0417 


0.0183 


-9.13 


• 10 


-3 


13.6 keV 


44.9 keV 


<j 7 : 263 MeV S/j : 3 


0.3444 


0.3450 


0.3752 


1.21 


• 10 


-2 


236 keV 


167 keV 


/i c (2P)^r?c7 

<j 7 : 821 MeV S/j : 3 


0.0646 


0.0631 


0.0489 


-5.44 


• 10 


-3 


161 keV 


142 keV 


X& 2 (2P) -> T(25)7 

<j 7 : 243 MeV 


—0.1186 


—0.1186 


—0.1221 


-1.38 


• 10 


-3 


18.5 keV 


16.4 keV 

Br: 16.2 ±2.4% 


X6i(2P) -> T(25)7 

q 1 : 229 MeV 


—0.1240 


—0.1241 


—0.1272 


— 1.23 


• 10 


-3 


16.8 keV 


15.1 keV 

Br: 21 ± 4 % 


Xm(2P) -> T(25)7 

q 1 : 207 MeV 


-0.1305 


—0.1306 


-0.1331 


—1.01 


• 10 


_3 


13.5 keV 


12.3 keV 

Br: 4.6 ±2.1% 


Xb2 (2P)^T 7 
q-y : 777 MeV 


-0.0176 


-0.0177 


-0.0156 


8.44 


• 10 


-4 


10.7 keV 


11.4 keV 

Br: 7.1 ± 1.0% 


Xbl(2P)^T 7 
q-y : 764 MeV 


-0.0155 


-0.0156 


-0.0132 


9.28 


• 10 


-4 


7.32 keV 


8.40 keV 

Br: 8.5 ±1.3% 


Yb0 (2P)^T 7 

q-y : 743 MeV 


-0.0123 


-0.0124 


-0.0098 


1.01 


• 10 


-3 


3.70 keV 


4.93 keV 

Br: 0.9 ±0.6% 


h b (2P)^r lb (2S) 1 

<? 7 : 257 MeV S fi : 3 


-0.1120 


-0.1120 


-0.1155 


-1.42 


• 10 


-3 


19.5 keV 


17.2 keV 


/i b (2P) -^7 

q-y : 821 MeV S/j : 3 


-0.0179 


-0.0179 


-0.0160 


7.43 


• 10 


-3 


13.4 keV 


13.5 keV 
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Table 9: The El dominated transitons between low-lying states in the B^ (c6) system. The labeling of 
the columns is as for Tables |]and[|. Note that spin-triplet states are indicated by "stars" in their labels. 
The q 1 values have been obtained as for Table ||. Further El transitions can be found in Table |l} 



Transition 






Mo [fin] 






M-2 


[fm] 




Width 




IA 


DYN 


El 








El 


DYN 


B* 2 -> B* 7 

q-t : 397 MeV 





1361 


0.1386 





1453 


2.43 


• 10 


-3 


i on VoV 
izu Ke v 


QQ Q LoV 
yo.ej Kt? V 


<j 7 : 379 MeV 





1382 


0.1408 





1469 


2.16 


• 10 


-3 


i n7 VoV 
± u ( Ke v 


OAR UV 


BSo - B* 7 

<j 7 : 352 MeV 





1403 


0.1429 





1480 


1.81 


• 10 


-3 


86.6 keV 


70.4 keV 


B* c (2S)->B* col 

<j 7 : 184 MeV 


-0 


1578 


-0.1578 


-0 


1610 


-1.28 


• 10 


-3 


4.71 keV 


4.22 keV 


B*(2S)^B* cll 

<? 7 : 157 MeV 


-0 


1720 


-0.1720 


-0 


1746 


-1.01 


• 10 


-3 


10.2 keV 


9.35 keV 


B*{2S) ~> B* 2 7 

<j 7 : 139 MeV 


-0 


1842 


-0.1843 


-0 


1864 


-0.85 


• 10 


-3 


13.2 keV 


12.3 keV 


B cl ^B cl 

<j 7 : 431 MeV S fi : 3 





1267 


0.1283 





1353 


2.50 


• 10 


-3 


135 keV 


103 keV 


B C {2S) -» S c i7 
<j 7 : 133 MeV 5/< : 3 


-0 


1834 


-0.1835 


-0 


1854 


-0.73 


• 10 


-3 


20.7 keV 


19.3 keV 



Table 10: The Ml transitions between low-lying S*-wave states in the Bf [cb) system. The matrix 
elements of the spin-flip operator associated with the OGE interaction are given in the column labeled 
"Oge". The 3S states have not been considered as they lie above the BD fragmentation threshold. 



Transition 


Matrix element [10 3 fm] 




Width 






NRIA 


RIA 


Conf 


Oge 


NRIA 


RIA 


RIA+Exch 


B*^B cl 

g 7 : 53 MeV 


18.51 


14.96 


-3.643 


3.969 


50.0 eV 


32.6 eV 


34.0 eV 


g 7 : 576 MeV 


1.015 


-1.437 


1.342 


1.182 


179 eV 


360 eV 


206 eV 


B*(2S)^B C (2S) 7 

<j 7 : 22 MeV 


18.49 


14.80 


-7.666 


2.480 


3.61 eV 


2.31 eV 


0.98 eV 


b c (2S)-b; 7 

q 7 : 507 MeV 


-1.067 


-3.089 


1.924 


8.351 


411 eV 


3.44 keV 


39.5 eV 
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Table 11: The El transitions from the 2P states in the B^r {cb) system. The labeling of the states and 
the photon momenta g 7 for the transitions are as for Tables |j and 



Transition 




M [fin] 




M 2 


[fm] 




Width 






IA 


DYN 


El 








El 


DYN 


B* c2 (2P)^ B* C {2S) 1 

<j 7 : 222 MeV 


0.2151 


0.2166 


0.2265 


3.94 


• 10" 


-3 


49.3 keV 


41.7 


keV 


B* cl (2P)^ B*(2S) 1 
g-y : 212 MeV 


0.2263 


0.2279 


0.2370 


3.59 


• 10" 


-3 


46.5 keV 


39.9 


keV 


5 t *o(2P)^ B* C (2S) 1 

q 1 : 194 MeV 


0.2381 


0.2398 


0.2474 


3.04 


• 10" 


-3 


39.0 keV 


34.2 


keV 


B* c2 {2P)^B* cl 

<j 7 : 733 MeV 


0.0320 


0.0333 


0.0264 


-2.97 


• 10" 


-3 


27.0 keV 


32.9 


keV 


B* cl (2P)^B* cl 

<j 7 : 723 MeV 


0.0270 


0.0283 


0.0204 


-3.34 


• 10" 


-3 


15.5 keV 


23.2 


keV 


B5,(2P)-B* 7 

<? 7 : 707 MeV 


0.0207 


0.0220 


0.0131 


-3.69 


• io- 


-3 


5.95 keV 


13.4 


keV 


B cl (2P)^B c {2S) 1 
q-y : 234 MeV S fi : 3 


0.2065 


0.2076 


0.2179 


4.08 


• io- 


-3 


53.3 keV 


44.5 


keV 


B cl {2P)^B cl 

<j 7 : 771 MeV S fi : 3 


0.0316 


0.0323 


0.0259 


-2.75 


• 10" 


-3 


30.6 keV 


35.7 


keV 



Table 12: The El transitions from the lightest spin-triplet D-wave states in the B^r (cb) system. The 
labeling of the states, the photon momenta g 7 and statistical factors Sfi are as for Table 



Transition 






Mo [fin] 




Width 








IA 


DYN 


El 


El 


DYN 


<j 7 : 235 MeV S fi : 


18/25 


0.2270 


0.2285 


0.2352 


75.7 keV 


65.4 


keV 


3 D 2 ^B* c2l 

q-y : 232 MeV S fi : 


9/50 


0.2279 


0.2294 


0.2358 


18.3 keV 


15.9 


keV 


3 D 2 B* cl 7 

qi : 250 MeV S fi : 


9/10 


0.2181 


0.2197 


0.2267 


63.9 keV 


54.7 


keV 


3 D t ^B* 2l 

<j 7 : 224 MeV Sfi : 


1/50 


0.2288 


0.2303 


0.2362 


1.84 keV 


1.61 


keV 


3 D, B* cl 7 

<j 7 : 243 MeV S fi : 


1/2 


0.2193 


0.2209 


0.2274 


32.5 keV 


28.0 


keV 


*D l ^B* c0l 

g 7 : 270 MeV S fi : 


2 


0.2079 


0.2096 


0.2170 


54.4 keV 


45.9 


keV 
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a c 


R* f 9 


R* f "i 


NRIA 


0.4810 n N 


0.4810 ^jv 


0.4810 fi N 


RIA 


0.4135 fx N 


0.4075 ^jv 


0.3973 ^jv 


Conf 


-0.0166 n N 


-0.0342 fi N 


-0.0487 fi N 


Oge 


0.0292 n N 


0.0186 ^jv 


0.0146 ^jv 


Total 


0.426 fi N 


0.392 n N 


0.363 //AT 



Table 13: The magnetic moments of the 
5- wave cb states in units of nuclear mag- 
netons. The RIA magnetic moment is 
given by eq. (fl9|), and the exchange cur- 
rent contributions from the scalar con- 
fining and OGE interactions by eqs. ( [sof ) 
and (jsi]), respectively. Note that the 
NRIA results are equivalent to those 
given by the static quark model. 



5 Discussion 



It is instructive to compare the numerical results of the previous section both with experiment and with 
those of other theoretical calculations, as there are several issues that are not readily apparent by casual 
inspection of the data presented in Tables ^ through [H^]. A direct comparison with the experimental 
averages given in ref. JTo[ ] may be misleading, as the branching fractions for various transitions are 
typically better known than the total width of the decaying state. Also, as recently published results j[lj 



for XcJ ~ * J 1^1 indicate that these El widths are significantly larger than previously thought 10 , then 
a review of the model predictions is called for. 



5.1 Ml transitions in the cc system 

• J/ip -^r\d 

The observed width of 1.14 ± 0.39 keV for the Ml transition J/ip — * rj c 7 has been difficult to explain 
theoretically since previous calculations typically overestimate it by a factor ~ 3. A possible 
solution for this overprediction, already hinted at in ref. ||, is presented in Table [| where the 
exchange current contribution from the scalar confining interaction brings the width down to the 
desired level. However, as shown in ref. juj, expansion of the RIA spin-flip operator to order v 2 /c 2 
overestimates the correction to the static quark model result, and in that case the usefulness of 
the exchange current contribution is not apparent. The importance of negative energy components 
for the transition J/ip — > rj c 7 has also been established within the instantaneous approximation 
to the Bethe-Salpeter equation in ref. jl2| and the Schrodinger approach in ref. |l4| , where widths 
close to that given in Table || were obtained for a scalar confining interaction. If the whole QQ 
potential had effective vector coupling structure, then no exchange current contribution would arise, 
as a vector interaction contributes a spin-flip operator only if the quark and antiquark masses are 
unequal, and agreement with experiment would thus be excluded. Other possible solutions include 
the introduction of a large anomalous magnetic moment for the charm quark || , but this possibility 
has apparently not been substantiated. 

This nonrelativistically forbidden Ml transition has also proved theoretically challenging, since the 
(near) orthogonality of the quarkonium wavefunctions leads to model-dependent results. In the 
recent calculation by ref. p2f, where good agreement with experiment was found for J/ip — > ?y c 7, 
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the width for ip' — > r] c 7 was however overpredicted by almost an order of magnitude. The present 
calculation gives a width of ~ 1.1 keV for that transition, which is close to the empirical value 
0.84 ± 0.24 keV [[ll]. This result is, however, not robust. That such a favorable result is obtained 
depends delicately on several factors, such as the employment of ip ' and rj c wavefunctions that 
include the spin-spin interaction in the S-w&ve, and the choice of approximation for the Ml matrix 
element. 

The amplitude ( |l5|) has the advantage of allowing the use of a realistic photon momentum in the 
expression (^3|) for the Ml width. It is useful to note that this treatment yields the same spin-flip 
operators than had the rigorous Ml approximation been used, as in the calculation of the exchange 
magnetic moment operators in refs. ^2f . Furthermore, the Ml approximation has been taken 
to affect the entire factor in brackets in eq. (|l5|). If the exponentials were separated from the 
current operators in eq. jl5|), then the width for ip ' — > r\ c 7 would be overpredicted by a factor ~ 4. 
However, if spin-averaged wavefunctions were employed, as in rcf. |T^| , then the conclusion would be 
exactly the opposite; In that case the present treatment would lead to unfavorable results. As seen 
from Table ^|, the exchange current operator associated with the scalar confining interaction gives 
the main contribution to the width for ip' — > r\ c 7 within this calculation. The present treatment of 
the Ml approximation may be regarded as consistent since it leads to the correct spin-flip operators 
and simultaneously allows the recoil of the r\ c to be taken into account. Further modifications to 
the Ml width for ip' — > rj c j may result from a relativistic treatment of the exchange current 
operator ( p3| ) and consideration of possible D-wave admixture in the ip' state due to the OGE 
tensor interaction. 

Ml transitions in the bb and cb systems 

T — > r] b 7 

A reasonable prediction for the width of the Ml transition T — ► rjb 7 can be obtained as soon as the 
mass of the rjb state is constrained by experiment, as the contributions from two-quark operators 
are suppressed by the large mass of the b quark. The unceratinty in the constituent mass (4885 
MeV) of the b quark is also nowadays rather small H [L8|. As realistic models of the spin-spin 
splittings for S- wave quarkonia give an rjb mass around 9400 MeV, then the width for T — > rjb 7 is 
likely to be less than 10 eV, as given in Table |[ However, the reliability of such predictions cannot 
be tested until the 775 state is discovered empirically. 

At this time, the photon momenta involved in calculations of Ml widths in the cb system have to 
be extracted from model calculations of the cb mass spectrum. As shown in ref. || , the mass of the 
spin triplet B* state is rather well constrained, and is expected to be about 6350 MeV. However, 
there is in general no such agreement for the magnitude of the B* — B c splitting, which determines 
the photon momentum of the B* — > B c 7 transition. Realistic models for the spin-spin interaction 
in the S'-wave appear to favor a small splitting of ~ 40 MeV. In spite of this, the computed width 
for B* — > £> c 7 of 34 eV given in Table [l(] compares well with the ~ 29 eV predicted by ref. || , 
even though exchange current contributions were not considered in that work. It is noteworthy 
that the exchange current contribution from the OGE interaction is dominant for cb and leads to a 
slight increase of the width for B* — > B c 7. Because of this cancellation, the nonrelativistic width 
is quite close to the net result. 
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Predictions have also been given, in Tables |^ and [Hj, for Ml transitions between QQ states that 
lie below the respective fragmentation thresholds. Notable among these is tp' — 77^7, which is similar 
to J/i/j — > r/ c j. However, the suppression due to the scalar confining interaction is stronger as the 2S 
wavefunctions have a longer range. Also, as recent experimental results indicate that the mass of the rj' c 
is larger than previously thought , then the phase space available for ip ' — > rf c 7 is also suppressed. 
Consequently the predicted width is also smaller than the values suggested by previous work . As for 
the transition ip ' — > rj c 7, the width for r)' c — > J/ip 7 is also sensitive to the particulars of the model because 
of cancellations in the matrix element and the large photon momentum. The results in Table |^ suggest 
that the width for this transition should be around 2 keV. As the experimental situation concerning the 
v[ c continues to improve, then the width for rj c — * J/ip 7 may possibly be measured in the near future. 

In the case of the bb system, the number of measurable Ml transitions is larger since the 35* states of 
bottomonium lie below the threshold for BB fragmentation. The results of Table || indicate that most 
Ml transitions in the bb system are difficult to calculate accurately, and therefore provide an instructive 
test for model predictions. In particular, the widths for transitions that do not change the principal 
quantum number of the quarkonium state are predicted to be highly suppressed, whereas the widths for 
transitions from excited r)b states to the T ground state are predicted to have larger widths of about 
100 eV. 



5.3 El transitions in the cc system 

While the Ml transitions in heavy quarkonia are sensitive to the Lorentz structure of the QQ interaction, 
the El transitions have been shown here to receive only small contributions from the exchange current 
operators of Fig. ^. They typically increase the value of the matrix element so that its value in the 
dynamical model is between those in the impulse and El approximations. On the other hand, the 
matrix element of the dipole operator is sensitive to the shape of the QQ wavefunctions, and thus a 
realistic description of the El widths requires that the hyperfine components of the QQ interaction are 
not treated as first order perturbations. This conclusion is in line with that reached in ref. fl8|| , which 
employed the nonsingular QQ potential model of ref. |32| . Comparisons between the present results, those 
of previous calculations and experiment are given for several El transitions in cc and bb in Tables [lj 
and [l| 



Table 14: Comparison of the results for El transitions in the charmonium (cc) system with the predictions 
of other models, for a scalar confining interaction. All widths are given in keV. The experimental widths 
have been extracted from ref. lO]. 





MB (ref. §) 


GS (ref. § H) 


MR (ref. §) 


This Work 


Exp (ref. 0) 


Xc2 -> 


J/ipl 


347 


413 


609 


343 


389 ± 60 keV 


Xcl -> 


J/ipl 


270 


340 


460 


276 


290 ± 60 keV 


XcO -> 


J/i>i 


128 


162 


225 


144 


165 ± 40 keV 




Xc2l 


27 


26 


41 


33 


20.4 ± 4.0 keV 




Xci 7 


31 


28 


48 


39 


25.2 ± 4.5 keV 




XcO 7 


19 


18 


37 


33 


26.1 ±4.5 keV 


h c — 


VcJ 


483 


630 




370 
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• XcJ -> J/lpl 

Predictions for the El transitions from the spin-triplet P-wave states in charmonium should be 
robust, as the quarkonium wavefunctions involved do not contain any nodes. The systematic over- 
prediction of these widths by ~ 20%, as indicated by the empirical data of ref. pTOj ] , has therefore 
been puzzling. However, in view of the most recent experimental data p] , this issue is apparently 
resolved. As seen from Table models which consider the recoil of the J/tp reproduce the empir- 
ical results remarkably well. In contrast, the rigorous El approximation leads to overpredictions 
of ~ 50%. 

• i> ' -> XcJ 7 

The El transitions from the ip' state have, in general, been difficult to predict, since the El 
approximation has typically overestimated the empirical widths by at least a factor ~ 2. The 
latest experimental data for the branching fractions for 7 decay and total width of the ip' suggests 
that the widths for ip' — > \ c j 7 should be around 25 keV, while the El approximation typically 
yields widths in excess of 40 keV. The results for ip ' — > XcJ 7 obtained by ref. jl^] within the 
framework of the instantaneous approximation to the Bethe-Salpeter equation, are also of the order 
~ 40 keV. Recoil effects cannot account for this overprediction, as they are much weaker than for 
the XcJ — * J/ip'y transitions. However, inspection of Table [l4| reveals that the predicted relative 
widths are also not in agreement with experiment, although the uncertainties are considerable. 
This suggests a sensitivity to small changes in the QQ wavefunctions, which is indeed revealed by 
inspection of the matrix elements in Table ^. The present empirical data suggests that the width 
for ip ' — > \c2 7 should be the smallest. Of the models presented in Table [u], the present one comes 
closest to this result. It should also be noted that significant reductions of the El widths have been 
obtained in ref. [[mJ by consideration of closed cq — qc fragmentation channels. 

The predicted widths for El transitions from the ip(3S) state have also been given in Table f|. The 
results suggest that the widths for transitions to the 2P states should be comparable to those for the 
ip ' — » XcJ 7 transitions. On the other hand, the widths for the ip{3S) — > XcJ 7 transitions are predicted to 
be smaller by a factor ~ 4. The empirical detection of any of these transitions will probably be challenging 
since the tp(3S) state decays mainly through DD fragmentation. However, the photon produced in the 
he ~* WcY transition will probably be detected in the near future, as that state lies well below the DD 
fragmentation threshold. The results presented in Table |IJ suggest that the El width for that transition 
is the largest in the cc system. However, the uncertainty in that prediction is somewhat larger because 
of recoil effects and the hitherto uncertain mass of the h c state. Of particular interest are also the El 
transitions from the Z D\ state given in Table [| as that state probably corresponds to the empirical 
i/>(3770) resonance. The calculated widths suggest that the transitions to the Xci and XcO states should 
be detectable by experiment, whereas that to the Xci state is highly suppressed by the statistical factor 
Sfi. 



5.4 El transitions in the bb system 

The number of measurable El transitions in the bb system is larger than in the cc system, as the T(3S') 
state lies below the threshold for BB fragmentation. While the branching fractions for several of the 
transitions presented in Table 15 have been measured, only very few of the total widths are known. This 
is due to the narrowness of the bb states, which makes a direct determination of their widths difficult. It 
is therefore instructive to compare the predicted El widths with those of other models, as well as with 
experiment. A review of the most important El transitions in Table hq is given below. 
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Table 15: Comparison of the predictions for El transitions in the bb system with those of other models 
that use a scalar confining interaction. All widths are given in keV. The experimental widths have been 
extracted from ref. p0[. Note that only the branching fractions are known for most of the transitions. 





/"i / p 111 ill \ 

GS (ref. g) 


GZ (ref. |8|) 


This Work 


Exp (ref. 


10) 


Xb2 " 


-> T 7 


33.0 


33.8 


36.0 


22 ± 3% 


Xbl ' 


-> T 7 


29.8 


30.4 


32.5 


35 ±89 




XbO - 


-T 7 


25.7 


25.3 


26.6 


< 6 % 




Jf 


> XbO 7 


0.73 


0.76 


1.01 


1.7 ±0.5 keV 


Jf 


> Xbl 7 


1.62 


1.37 


1.80 


3.0 ± 0.7 keV 




► Xb2 7 


1.84 


1.45 


2.03 


3.1±0.7keV 


X' b2 - 


-»• T'7 


10 

iz.y 


1 £ 9 

10. z 


ID. 4 


16.4 ±2.4 % 


Xbi - 


-+ T'7 


11.9 


14.7 


15.1 


21 ± 4 % 


Xbo ' 


-> T'7 


10.6 


12.3 


12.3 


4.6 ±2.1 % 


Xb2 " 


-^T 7 


18.2 


10.4 


11.4 


7.1 ± 1.0 % 


X'bi - 


— > T 7 


11.8 


7.51 


8.40 


8.5 ± 1.3 % 


x' b o - 


-T7 


6.50 


3.57 


4.93 


0.9 ±0.6 % 


Y« _ 


+ Xw 7 


0.114 


0.029 


0.15 






"£ll _ 


+ Xw 7 


0.003 


0.095 


0.11 






y// _ 


+ Xb2 7 


0.194 


0.248 


0.04 






^£11 


*Xbo7 


1.09 


1.30 


1.14 


1.4 ±0.3 keV 




* Xbi 7 


2.15 


2.34 


2.12 


3.0 ± 0.5 keV 


y// _ 


^X 62 7 


2.29 


2.71 


2.50 


3.0 ± 0.6 keV 



• Xbj -> T 7 

The calculated widths for Xbj — > T7 agree rather well with those of the previous calculations 
presented in Table |l5|, although they are generally somewhat larger. The results of the present 
calculation suggest that the total width of the Xb2 state should be 164 ±22 keV and that of the Xbi 
93 ± 22 keV. Similarly, the total width of the Xbo is probably larger than <~ 440 keV. This situation 
is similar to that observed for cc |Tcj , where the Xc2 is wider than the Xci by about a factor ~ 2. 

• T' -> Xbj 7 

The largest uncertainty concerning the T(2S f ) — > XbJl is the measurement of the total width 
of the T(2S) state. Originally given as ~ 27 keV, it now stands at 44 ± 7 keV 0. Thus the 
model predictions in Table |lj, which originally fitted the experimental data well, no longer do so 
satisfactorily. It is therefore very difficult to judge the quality of any given prediction until the 
experimental situation is stabilized. Still, it is noteworthy that the present calculation does give 
slightly better agreement with experiment than the previous models in Table [l^. 
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• Transitions from the x'bj states 

The El transitions from the Xbj(2P) states in bottomonium provide a useful test for theoretical 
models since experimental data now exists on all six branching fractions |lC| ] , even though the total 
widths of the Xbj{^P) states are not known. The experimental results indicate that the widths for 
transitions to the T should be about one half of those for transitions to the T(25), even though 
much more phase space is available for the former. Indeed, it can be seen from Table [Hi] that 
the model of rcf. |3|], where spin-averaged wavefunctions were employed, does not compare very 
well with the experimental branching fractions even though the hyperfine splittings of the Xbj(2P) 
states are quite small. A realistic description of these transitions requires that the hyperfine effects 
are accounted for by the QQ wavefunctions, in which case the computed El widths agree much 
better with experiment. 

As the calculated widths for Xbj(2P) — » T(25)7 are almost the same in ref. jy| and the present 
model, then it is possible to obtain realistic estimates for the total widths of the Xbj(2P) states 
from the measured branching fractions for the El transitions Xbj(2P) — > T(25)7. The width of 
the Xb2^P) state is then predicted to be 100±15 keV, while that of the xw(2F) is 72±14keV. The 
Xbo(2P) state appears to be significantly broader, but because of the large errors in the reported 
El branching fractions, only a rough estimate of 267 ± 140 keV is possible. 

• Transitions from the T" state 

As the reported total width of the T(35) state |lCj, 26.3 ± 3.5 keV, is better known than that 
of the T(2S) state, then it is expected that systematic uncertainties in the reported experimental 
results for T(3S) — > Xbj(%P) 7 should be smaller than for the analogous T(25) — ► Xb.J 7 transitions. 
By inspection of Table [l^, it can be seen that the T(35) — > Xbj{^P)l transitions are generally 
rather well described by a number of models, although the calculation of ref. where spin- 
averaged wavefunctions were employed, appears to underpredict the empirical widths. Also, the 
results of ref. appear to compare slightly more favorably with experiment than does the present 
calculation. 

On the other hand, the situation concerning the T(3<S) — ► XbJl transitions remains unsettled 
because of a strong cancellation in the El matrix element. As all of the models presented in Table |lj 
predict different widths for the T(3S I ) — > Xbj 7 transitions, then only experimental determination 
of the widths can settle the question. However, this may turn out to be a formidable task since all 
models predict widths that are an order of magnitude smaller than that of any previously measured 
El transition in the bb system. It is seen by inspection of Tabled that the calculated widths in the 
El approximation are similar to those of ref. || in that the width for T(3S I ) — * Xbl 7 is vanishingly 
small. However, the dynamical model predicts that the width for T(35) — ► Xbo 7 should be the 
largest and that for T(35) Xb2j the smallest. It is encouraging that the same pattern is also 
predicted in Table || for the analogous transitions in the cc system, where the widths are much 
larger relative to the other El transitions. It is also noteworthy that the T(3S) — * Xbj 7 transitions 
obtain an appreciable contribution from the matrix element A^2 in Table pi 

In addition to the transitions considered in Table [l^, the El decays of the spin-singlet ht and 77^ 
states and the lightest spin-triplet D-wave states have also been calculated in Tables ||, fj] and ||. No 
empirical data exists as yet on any of these states. The pattern of El widths for these states is predicted 
to be similar to that for the analogous states in the cc system, although the widths are much smaller for 
bb because of the narrower wavefunctions involved. 
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5.5 El transitions in the cb system 



The determination of photon momenta for El 
transitions in the bottom-charm Bf mesons 
has to rely completely on model predictions 
for the masses of the cb states. However, 
the uncertainty introduced by this is small, 
as model predictions for the major level 
splittings agree with each other to a large 
extent g . The results in Table |l| reveal that 
the predictions of the present model are sim- 
ilar to those obtained by ref. H, although 
significant differences exist for transitions 
such as B*(2S) -> B* j and B* 2 (2P) -> 
B*(2S) r y, where the widths are sensitive to 
the effects of the hyperfine components of the 
QQ interaction. It is noteworthy that while 
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the predicted widths for the B* cl — > B* c 
transitions agree rather well with those from 
ref. Q , there is a significant disagreement for 
B c i — > -B c 7. When the somewhat different 
photon momenta are accounted for, this dis- 
agreement amounts to about a factor ~ 3. 
An issue not considered in this paper is the 
mixing of the L = 1 states with J = 1 , which 
is due to the antisymmetric spin-orbit inter- 
action that was not included in the Hamilto- 
This mixing, which was considered 
has the effect of allowing " spin- flip" 
El transitions of the type B* cl — > B c ^. How- 
ever, the widths for such "forbidden" tran- 
sitions were found in ref. Q to be typically 
suppressed by a factor ~ 100 relative to the 
"allowed" ones considered in this work. 




Table 16: Comparison of the results for El transitions 
in the B^r (cb, be) system with the predictions of ref. H . 
All widths are given in keV. In the notation B*j, where 
J is the total angular momentum, states without stars 
in their labels are spin singlets in the LS coupling 
scheme. 
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As the magnetic moment operators of a QQ system have been derived in this work, then the magnetic 
moment of the S*-wave spin-triplet B* states have been calculated as an interesting by-product. It is 
seen from Table that the magnetic moments also receive significant corrections from the two-quark 
operators considered in this paper. This is not surprising since the situation is similar for the Ml decays 
of the QQ systems. Inspection of Table [l^ reveals that the net relativistic decrease of the B* magnetic 
moment amounts to about 15 %, and that the exchange current contributions actually increase the net 
magnetic moment slightly. This situation was also noted for the Ml transition B* — » _B C 7. However, 
because of the short-range nature of the OGE interaction, its contribution quickly becomes subdominant 
for the higher 5-wave states. The predicted magnetic moments of the B* mesons are in line with the 
work of ref. p3| on the magnetic moments of the baryons. 
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5.6 Other frameworks 



The masses, radiative transitions and other properties of the heavy quarkonia have also been investi- 
gated within the frameworks of QCD sum rules ||, fjj, Heavy Quark Effective Theory (HQET) || and 
Non-Relativistic QCD (NRQCD) ||. It is of particular interest to compare the results of the present 
calculation with those of ref. |^|, where the El transitions between cc and bb states were considered. The 
results of ref. Q are presented in terms of an "effective" matrix element 5, from which the dependence 
on the angular momenta of the QQ states and the momentum of the emitted photon has been extracted. 
The theoretical uncertainties are therefore concentrated into the parameter S. As that parameter is 
closely related to the radial matrix element Mo for El decay, a direct comparison between the present 
work and that of ref. JsJ is possible. Such a comparison is given in Table [l7]. In addition, the width for 
the El transitions h c — > r/ c "f was in ref. |s) obtained as ~ 450 keV, which in view of the somewhat larger 
photon momentum of ref. [0 compares reasonably well with the present value of 370 keV. 



Table 17: Comparison between the effec- 
tive matrix element 5 of ref. |^] and the 
corresponding values extracted from the 
present calculation. The exact definition 
of 6 can be found in ref. ||. Note that in 
the present model, the matrix element for 
an El transition also depends on the total 
angular momentum J, as spin-averaged 
wavefunctions are not employed. The re- 
ported values in this table therefore cor- 
respond to the spin average of the matrix 
elements given in Tables [| and ||. 
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+ XcJl 


0.209 ±0.012 


0.267 


Xc.J - 


-> J/-07 


0.197 ±0.012 


0.215 




+ XbJl 


0.109 ±0.009 
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J// 
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0.152 ±0.009 
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5.7 Conclusions 

The conclusion of this work concerning the exchange current operators and the associated two-quark El 
and Ml operators is that they play a major role in spin-flip Ml transitions, whereas they are insignificant 
for the El transitions. The reason for the smallness of the exchange charge contributions to the El widths 
is the large masses of the charm and bottom quarks, as those contributions are proportional to m~ 3 . 
However, as suggested in ref. [[35|, such operators may be significant when light (~ 400 MeV) constituent 
quarks are involved. 

A satisfactory description of the Ml widths of charmonium was achieved in this work, under the 
assumption that the effective confining interaction is purely scalar. This is in line with the results of 
ref. |E| , where the width for J/ip — > rj c 7 was found to be well described by an effective scalar confining 
interaction within the framework of the instantaneous approximation to the Bethe-Salpeter equation. 
This is reassuring since it has been shown in ref. |3(| that models (of equal mass quarkonia) which 
employ an effective linear vector confining interaction or a superposition of scalar and vector confining 
interactions with positive weights are inconsistent with the properties of QCD. 

It should be noted that the instanton induced interaction for Qq and QQ systems, as given by ref. p7| , 
may also contribute a significant two-quark exchange charge operator. For light constituent quarks, such 
contributions may be large whereas they are much smaller for charm quarks pp| |. A calculation of 
the associated spin-flip operators for Ml transitions could therefore provide useful and constraining 
information on the strength of the effective instanton induced interaction in mesons with heavy quarks. 
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